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Abstract 


Introduction 


Repeated  solutions  to  the  continuity 
and  vorticity  equations  are  frequently 
required  in  computations  of  three  dimen- 
sional flows  in  turbomachinerv  passaqes. 
When  the  two  equations  are  nonhomoqeneous 
previous  formulations  resulted  in  two 
second  order  differential  equations.  A 
new  approach  is  presented  here,  which  is 
applicable  in  a generalized  two  dimensional 
domain  or  axisymmetric  field.  It  is  based 
on  the  definition  of  a streamlike  function 
which  is  used  to  transform  these  nonhomo- 
geneous  first  order  partial  differential 
equations  to  a single  second  order  equation 
with  Dirichlet  boundary  conditions  over  the 
solid  boundaries.  Some  applications  are 
presented  to  show  how  this  new  approach  can 
be  used  to  save  computer  time  in  numerical 
flow  solutions. 


Nomenclature 


C,  to  Cc 

I D 

j 

s 

u 

v 

X 

r 

yr 

and  a ^ 


arbitrary  constants  in  Eqs. 

(11)  to  (16). 

parameter  equal  to  zero  for  two 
dimensional  flow,  and  to  one  for 
axisymmetric  flow 

source/sink  term  in  Eq.  (1) 
velocity  component  in  x-direction 
velocity  component  in  y-direction 

reference  value  of  the  x coor- 
dinate in  the  x^  formulation 

reference  value  of  the  y coor- 
dinate in  the  X2  formulation 

the  nonhomogeneous  terms  in  the 
second  order  equations  for  the 
streamlike  functions,  Eqs.  (5), 
(6),  (9)  and  (10) 


x.  and  x5  streamlike  functions  defined  in 
1 Eqs.  (3),  (4),  (7)  and  (8) 

w flow  vorticity,  Eq.  (2) 


The  simultaneous  solution  of  the  two 
first  order  partial  differential  equations 
representing  the  conservation  of  mass  and 
the  vorticity,  is  required  in  many  flow 
studies. The  authors  interest  in  this 
problem  is  connected  to  their  internal 
nonviscous  subsonic  flow  investigations ^ i 
in  the  various  turbomachinery  passages. 

The  need  for  the  solution  to  the  outlined 
problem  is  also  encountered  in  other 
diversified  flow  fields  such  as  external 
two  dimensional  transonic  flows^  and 
internal  viscous  flows.**'5  In  three 
dimensional  flow  fields,  the  equation  of 
conservation  of  mass  includes  a derivative 
of  the  third  velocity  component.  There- 
fore, except  for  two  dimensional  incom- 
pressible flow  or  for  irrotational  flow, 
the  two  equations  are  generally  nonhomo- 
geneous.  When  at  least  one  of  the  equations 
is  homogeneous,  the  traditional  formulation 
of  flow  problems  has  been  in  terms  of  a 
potential  function  or  a stream  function. 

In  the  absence  of  a source  term  in  the 
continuity  equation,  the  stream  function 
is  introduced  into  the  rotationality  equa- 
tion to  obtain  the  second  order  partial 
differential  equation  for  the  stream 
function.  On  the  other  hand,  the  irrota- 
tional flow  is  usually  studied  in  terms  of 
the  potential  function  which  is  used  into 
the  continuity  equation  to  obtain  the 
governing  second  order  partial  differential 
equation.  There  has  been  no  single  unified 
approach,  however,  for  the  mathematical 
formulation  of  the  problem  when  both  the 
continuity  and  the  rotationality  equations 
are  nonhomogeneous.  In  references  4,  5 
and  6,  two  different  mathematical  approaches 
were  used  for  the  solution  of  the  con- 
tinuity and  rotationality  equations.  Rubin 
and  Khosla^  formulated  the  problem  in  terms 
of  a stream  function  and  a potential 
function,  resulting  in  two  second  order 
partial  differential  equations  in  these 
functions.  The  problems  of  handling  the 
Neumann  boundary  conditions  for  the  poten- 
tial function  over  irregular  boundaries 
are  well  known. ? There  are  additional 
nroblems  involved  in  the  numerical  solution 
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using  finite  difference  methods,  when  the 
potential  equation  has  a nonhomogeneous 
term,  with  this  type  of  boundary  conditions 
over  all  the  boundaries  of  the  solution 
domain.  > 8 A different  approach  was 
followed  in  references  5 and  6.  Cross 
differentiation  was  used  to  obtain  two 
second  order  partial  differential  equations 
in  the  two  velocity  components.  In  this 
case,  additional  boundary  conditions  are 
required  which  are  obtained  from  the 
original  two  first  order  equations. 

The  authors  are  proposing  in  this 
paper  a new  unified  approach  for  the  mathe- 
matical formulation  of  the  nonhomogeneous 
continuity  and  rotationality  eouations. 

The  approach  is  based  on  the  definition 
of  a new  dependent  variable  which  will  be 
referred  to  as  "a  streamlike  function". 

The  introduction  of  this  new  dependent 
variable  transforms  the  two  first  order 
nonhomogeneous  partial  differential  equa- 
tions to  a single  second  order  equation  in 
the  new  variable.  This  formulation  is 
clearly  superior  therefore  to  those  used  in 
references  4,  5 and  6,  which  resulted  in 
two  simultaneous  second  order  equations. 

The  advantage  of  the  new  formulation  in 
terms  of  computer  time  savings  are  obvious, 
since  only  half  the  number  of  the  dif- 
ference equations  need  to  be  solved. 

Another  advantage  of  the  formulation  is 
that  the  boundary  conditions  for  the  stream- 
like function  are  of  the  Dirichlet  type 
over  the  solid  boundaries.  This  can 
particularly  be  very  helpful  in  the  numeri- 
cal solution  of  flow  problems  with  irregular 
boundaries,  using  finite  difference  methods. 
Some  applications  of  the  streamlike 
function  formulation  to  the  flow  in  turbo- 
machine passages  will  be  presented  to  show 
the  diversity  and  flexibility  of  this  new 
approach. 


follows : 


1 3xl 


u = -j  a^r  + -j  / xD  s<x'y> 

V J ' V J v 


Where  the  subscript  r refers  to  a chosen 
reference  value. 

When  equations  (3)  and  (4)  are  sub- 
stituted into  equation  (2),  we  obtain: 


3 Xi  3 x 


i A1 

t-  ~ x sir  - °i(x'y) 


° ^ (x , y)  = x3  u(x,y)+|—  / x3  S(x,y)  dx 


The  deviation  from  the  standard  definition 
of  the  stream  function  in  this  case  is 
in  the  x-velocitv  component  u,  given  by 
eauation  (3).  Another  streamlike  function, 
X2<  that  satisfies  eauation  (1)  automati- 
cally, could  also  be  defined  as  follows: 


3^-  + -j  f x3  S(x,y)  dy 


Mathematical  Formulation 

The  continuity  and  rotationality 
equations  can  generally  be  written  in  the 
following  form: 

( x3 u)  + |^  (x3v)  = xj  S<x,y)  (1 


When  this  definition  is  substituted  into 
eauation  (2) , the  following  equation  is 
obtained  for  the  streamlike  function  X2: 


32X2  32x2  j 3x2 

3X2  + 3y2  ‘ X 3X 


- o 2 (x,y) 


!y  <u)  " h.  (V)  = " “(x'y) 


Where  u,  v are  the  velocity  components  in 
the  x and  y directions  respectively,  S is 
the  source/sink  term,  and  ui  is  the  vorti- 
city  in  the  direction  normal  to  plane  x,v. 
The  value  of  j is  zero  for  two  dimensional 
flow  and  j is  one  for  the  axisymmetric 
case. 

The  Streamlike  Function 

A new  streamlike  function  is  intro- 
duced, such  that  the  continuity  equation 
is  automatically  satisfied.  The  velocity 
components  are  defined  in  terms  of  the 
streamlike  function  and  the  source  term  as 


o.(x,y)  = x3u(x,y)  - xJ  — / S(x,y)  dy 

2 3x  yr 

(10) 

In  this  case,  the  deviation  from  the 
traditional  stream  function  appears  in 
equation  (8),  which  defines  v,  the  velocity 
component  in  the  v direction. 

It  can  be  seen  from  eauations  (6)  and 
(10)  that  the  nonhomoaeneous  terms  oj  and 
o2  in  the  resulting  second  order  equations 
are  not  only  dependent  on  the  vorticity  w, 
but  are  also  dependent  on  the  source  term 
S in  the  original  continuity  equation,  and 
on  the  choice  of  the  reference  coordinate 

x or  y . 
r 1 r 


The  Boundary  Conditions 


Two  types  of  boundary  conditions  will 
be  discussed,  namely  that  involvinq  a 
specified  velocity  component  in  the  di- 
rection perpendicular  to  a boundary,  and 
the  other  involving  a specified  velocity 
component  in  the  direction  tangent  to  a 
boundary. 


Specified  Normal  Velocity: 

This  type  of  boundary  condition  is 
encountered  when  the  volume  flux  rate  is 
specified  over  a given  portion  of  the 
boundary.  Zero  normal  velocity  components 
are  usually  associated  with  stationary 
impermeable  solid  walls. 

i.  u = f(y)  on  x = constant,  can  be 
expressed  in  terms  of  x^  as 

. x 

X^  = / X3  f(y)  dy  - //  x-'stxjy)  dxdy  + C1 


and  can  be  expressed  in  terms  of  x2  as 
X 2 = / xj  f (y)  dy  + C2  (12) 


ii.  v = g(x)  on  y = constant,  can 
be  expressed  in  terms  of  Xj  as 


^ = / x-1  g ( x)  dx  + C., 


(13) 


and  can  be  expressed  in  terms  of  X2  as 


= - S x-’  g (x)  dx  + //  x-^S  (x,y)  dydx  + C4 

yr 

r (14) 


iii.  More  generally,  when  the  normal 
velocity  component  Vn  is  specified  over  a 
general  irregular  boundary  £,  the  boundary 
conditions,  in  terms  of  the  streamlike 
function,  x,,  can  be  expressed  as: 


V d£  - 
n 


x 

ff 

X 

r 


x-’sixjy)  dxdy  + C5 


(15) 


and  in  terms  of 
X2  as 

X 2 = / xj  Vn  d£ 


the  streamlike  function 
+ //  xJS(x,y)  dydx  + Cg 


where  Cj  to  Cg  in  equations  (11)  to  (16) 
represent  arbitrary  constants. 

It  is  clear  from  equations  (11) 
through  (16)  that,  Dirichlet  type  boundary 
conditions  for  the  streamlike  function, 
result  when  the  flow  velocity  normal  to 
the  boundary  is  specified.  As  expected, 
the  line  integral  of  the  volume  flux 


rates  contributes  to  the  variation  in  the 
streamlike  function  over  a boundary.  In 
addition,  the  area  integrals  of  the  source 
terms  S(x,y)  in  equations  (11),  (14),  (15) 
and  (16)  can  account  for  streamlike  function 
variation  over  a boundary,  through  which 
there  is  no  flux. 


Specified  Tangential  Velocity: 

In  this  case  the  boundary  conditions 
are  of  the  Neumann  type  when  expressed  in 
terms  of  the  streamlike  function. 


i.  v = q(y)  on  x = constant,  can  be 
expressed  in  terms  of  Xj  us 


q (y) 


(18) 


and  can 


3x 


be  expressed  in 

i y 

- x J g (y)  + ! 


yr 


terms  of  x2  as 
x^SU.y)  dy 


(19) 


ii.  u = f(x)  on  y = constant,  can  be 


expressed  in  terms  of  X^  as 

3xl  j X 1 

•r—  = xJ  f (x)  - / xJS(x,y)  dx 
y v 


(20) 


and  can  be  expressed  in  terms  of  x2  as 


3X 

9 V 


2 - xj  f(x) 


(21) 


This  type  of  boundary  conditions  will  not 
be  formulated  over  a general  boundary  shape, 
since  it  is  generally  associated  with  the 
boundary  conditions  at  infinity.  Fquations 
(18)  through  (21)  are  sufficient  for  this 
purpose . 


Numerical  Methods 


The  rest  of  this  paper  is  mainly 
intended  for  showing  how  the  new  streamlike 
function  formulation  can  be  used  in  con- 
junction with  already  existina  numerical 
methods  for  second  order  equations  to 
obtain  solutions  to  flow  problems.  A 
brief  review  of  the  available  numerical 
methods  for  solving  nonhomogeneous  second 
order  differential  equations  is  appropriate 
therefore  at  this  point.  Before  proceeding 
with  this  review  however,  we  will  describe 
very  breiflv  the  numerical  methods  which 
are  available  for  solving  the  two  first 
order  equations  directly  for  the  velocity 
components.  These  methods  are  based  on 
writing  the  finite  difference  form  of  these 
eauations  at  staqgered  arid  points  rather 
than  at  the  same  arid  point.  This  idea 
was  first  introduced  by  Gates  and  Von 
Rosenberg! 0 who  used  centered  difference 
schemes  for  expressing  the  first  order 
derivatives  of  the  flux  components  in 


3 


potential  flow.  They  developed  a direct 
numerical  method  and  various  implicit  line 
iterative  methods  for  solving  the  resulting 
sets  of  equations  in  the  flux  components. 
Using  different  grid  sizes,  they  evaluated 
the  accuracy  and  the  approximate  total 
number  of  arithmetic  operations  involved 
in  each  of  their  methods  of  solution  in  a 
given  regular  domain.  Martin  and  Lomax** 
employed  the  same  difference  scheme  over  a 
staggered  grid  in  their  solution  of  the 
nonhomogeneous  continuity  and  rotat ional itv 
equations.  In  their  procedure,  they  mani- 
pulated the  resulting  finite  difference 
equations,  to  obtain  a set  of  algebraic 
equations  in  only  one  of  the  velocity  com- 
ponents which  they  solved  using  the  cyclic 
reduction  method. 

Relaxation  methods  have  been  used  for 
a long  time  and  are  still  in  great  use  for 
the  solution  of  elliptic  equations.*"* 
Recently,  several  investigators  became 
interested,  however,  in  the  development  of 
fast  direct  methods  for  solving  the  cen- 
tered difference  form  of  Poisson's  equa- 
tion. The  first  two  of  these  investiga- 
tions which  were  reported  by  Buneman**  and 
Hockney* ^ use  the  cyclic  reduction  method, 
and  the  finite  Fourier  transform  method, 
respectively,  for  the  numerical  solution 
of  Poisson  equations  with  Dirichlet  bound- 
ary conditions.  Several  other  investiga- 
tors*2»13  have  contributed  since  then  to 
reducing  the  limitations  to  the  application 
of  these  two  methods.  The  improvements 
included  generalizing  Buneman's  cyclic 
reduction  method  to  Neumann  and  periodic 
boundary  conditions  on  regular  boundaries 
and  to  be  used  with  arbitrary  number  of 
grids  by  Sweet. *^  Hockney's  finite  Fourier 
transform  method  was  applied  to  problems 
with  Dirichlet  boundary  conditions  on  irreg- 
ular boundaries,  and  with  Neumann  boundary 
conditions  on  regular  boundaries  by  Buzbee 
et  al.*6  in  cases  involving  complicated 
boundary  conditions,  the  relaxation  methods 
have  remained  to  be  the  alternative  method 
of  solution.*7'*8  Although  these  methods 
are  generally  much  simpler  to  program,  it 
has  been  found  that  the  boundary  conditions 
have  to  be  carefully  handled  in  the  case 
of  nonhomogeneous  elliptic  equations  with 
Neumann  boundary  conditions  if  a convergent 
solution  is  to  be  obtained.2'8'9 


Applications 

Two  example  problems  are  presented  to 
illustrate  the  use  of  the  streamlike  func- 
tion in  the  solution  of  flow  problems  in 
turbomachinery  passages.  The  first  set  of 
results  simulate  the  secondary  flow  in  the 
planes  perpendicular  to  the  through  flow 
direction  for  the  nonviscous  rotational 
flow  in  a curved  duct.  The  numerical  pro- 
cedure for  determining  this  three  dimension- 
al flow  field  was  described  in  reference  1; 
here  we  only  present  the  results,  usinq 
the  present  formulation  for  a given  vorti- 
city  and  source  distribution.  The  case 
under  investigation  represents  a source 
distribution  with  linear  variation  in  the 


horizontal  direction,  with  a mean  value  of 
zero  and  maximum  and  minimum  values  of  +1 
and  -1  at  the  walls.  The  vortex  distribu- 
tion changes  linearly  in  the  vertical 
direction  with  a maximum  value  of  1.5  at 
the  bottom  wall  and  a minimum  value  of  0.5 
at  the  upper  wall.  The  results  for  this 
case  are  presented  in  Figures  1 and  2 for 
the  first  and  second  streamlike  function 
formulations,  x^  and  > respectively. 

When  xr  is  taken  outside  the  solution 
domain,  the  xi  contours  are  symmetric  about 
the  vertical  centerline  as  shown  in  Fig.  1 
and  do  not  change  whether  xr  is  taken  to 
the  right  or  to  the  left  hand  side.  With 
u = 0 over  the  boundaries  x = constant  and 
v = 0 over  the  boundaries  y = constant, 
equations  (11)  and  (13)  give  constant 
values  of  xi  over  all  the  solid  boundaries 
as  shown  in  Fiq.  1.  Furthermore,  the 
streamlike  lines,  xj  = constant,  constitute 
closed  contours  inside  the  solution  domain 
in  this  case.  In  the  second  formulation, 

X2  *s  also  constant  over  the  two  solid 
boundaries  x = constant  and  the  X2  contours 
are  symmetric  about  the  vertical  centerline. 
The  streamlike  function  X2  *s  not  neces- 
sarily constant  however  over  the  two 
boundaries,  y = constant.  When  yr  is 
chosen  outside  of  the  solution  domain,  X2 
remains  constant  over  the  wall  facing  the 
reference  line,  but  varies  along  the 
opposite  wall.  This  means  that  the  stream- 
like lines  will  intersect  that  boundary 
opposite  to  the  side  where  yr  is  chosen. 

This  is  seen  in  Fig.  2,  in  which  y,-  was 
chosen  on  the  lower  side.  Although  the 
streamlike  lines  are  dependent  on  the 
formulation  and  the  choice  of  the  reference 
line,  the  actual  flow  velocities  are  inde- 
pendent of  these  choices.  The  arrows  in 
Fig.  3 show  the  magnitude  and  direction  of 
the  resulting  secondary  flow.  While  the 
previous  results  were  for  a domain  with 
regular  boundary,  the  rest  of  the  appli- 
cations are  for  a flow  problem  with  an 
irregular  boundary. 

The  second  set  of  results  are  pre- 
sented for  the  nonviscous  incompressible 
irrotational  flow  in  the  cross-sectional 
planes  of  a radial  inflow  turbine  scroll. 

In  this  case,  the  vorticity  u in  the 
rotationality  equation  is  equal  to  zero 
everywhere  in  the  flow  field,  and  the 
source  term  S,  in  the  continutiy  equation 
is  dependent  on  the  through  flow  velocity 
profile  variation. ^ The  boundary  conditions 
for  this  case  consist  of  specified  normal 
velocity  component,  over  all  the  boundaries 
of  the  solution  domain  shown  in  Fig.  4. 

The  velocities  normal  to  the  solid  boundary 
ABCD  and  to  the  axis  of  symmetry  EA  are 
equal  to  zero,  and  the  velocity  normal  to 
the  scroll  exit  DE  is  uniform  and  different 
from  zero.  In  reference  2 this  problem 
was  formulated  in  the  traditional  way, 
using  the  potential  function,  which  re- 
sulted in  a Poisson  equation  with  Neumann 
boundary  conditions.  Because  of  the 
irregular  boundary  shape,  the  numerical 
solution  was  obtained  using  relaxation 
methods.  A large  number  of  iterations  were 
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required  before  the  solution  converged  in 
spite  of  the  special  care  in  handling  the 
boundary  conditions.  The  same  problem  is 
reformulated  here  in  terms  of  the  stream- 
like  functions  with  the  results  shown  in 
Figs.  5 through  8 for  the  case  with  the 
uniform  source  distribution.-  In  all 
cases,  the  value  of  the  streamlike  function 
was  taken  equal  to  zero  at  the  corner 
point  A.  Both  Fios.  5 and  6 represent  the 
results  for  the  first  streamlike  function 
formulation,  \j.  The  difference  between 
the  resulting  streamlike  lines  in  the  two 
figures  is  due  to  the  different  choices  of 
the  reference  x-coord inate , xr.  In  the 
first  case,  xf  was  taken  external  to  the 
solution  domain  to  the  right  hand  side,  and 
in  the  second  case,  xr  was  taken  external 
to  the  solution  domain  to  the  left  hand 
side.  According  to  equation  (15) , when 
xf  is  on  the  right  hand  side,  the  stream- 
1 ike  function  xj  remains  constant  over  the 
part  BC  of  the  curved  boundary  as  shown 
in  Fig.  5,  and  when  xr  is  to  the  right, 
remains  constant  over  the  part  AB  of  the 
curved  boundary  as  shown  in  Fig.  6.  In 
both  cases,  the  streamlike  function  xi 
remains  constant  over  the  straight  boundary 
portions  CD  and  AE  according  to  equation 
(13)  since  the  flow  velocity,  v,  is  equal 
to  zero. 

Figures  7 and  8 represent  the  results 
obtained  using  the  second  streamlike 
function  formulation,  with  vr  taken 

above  the  solution  domain  in  the  first 
case  and  below  the  solution  domain  in  the 
second  case.  It  is  interesting  to  notice 
that  when  yr  is  taken  below  AE , the  right 
hand  side  of  equation  (10)  becomes  identi- 
cally equal  to  zero.  As  in  the  previous 
formulation,  the  streamlike  function,  \j, 
remains  constant  over  the  solid  boundaries 
to  the  side  on  which  yr  is  chosen  outside 
the  domain. 


Discussion 

It  has  been  shown  through  the  results 
presented  that  the  choice  of  the  type  of 
streamlike  function  formulation  and  the 
reference  line  alter  not  only  the  boundary 
conditions,  but  also  the  nonhomogeneous 
term  in  the  resulting  differential  equation. 
This  choice  can  be  used  to  obtain  the 
boundary  condition  which  is  easiest  to 
handle.  This  can  usually  be  accomplished 
by  placing  the  reference  line  on  the  side 
of  the  most  irregular  solid  boundary  to 
have  constant  value  of  the  streamlike 
function  over  it.  The  numerical  solution 
itself  can  be  sensitive  to  these  choices, 
as  the  authors  found  in  the  case  of  the 
scroll  where  the  flow  velocities  are  very 
small  in  the  part  opposite  the  exit  neck. 
Taking  the  reference  line  to  the  right, 
with  the  xj  formulation,  was  found  to 
result  in  the  closest  spaced  stream- 
like lines  in  this  region.  In  all  the 
results  presented  here,  the  reference 
lines  were  generally  chosen  external  to 
the  solution  domain  for  the  purpose  of 
demonstration.  The  formulation  itself 


does  not  place  any  restriction,  however, 
on  placing  the  reference  line  inside  the 
solution  domain.  The  choices  for  the 
location  of  the  reference  line  in  any 
situation  are  infinite  and  should  be  deter- 
mined by  the  user  depending  on  the  type  of 
his  problem. 
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Streamlike  Function,  \j,  Contours. 


Fig.  2.  Streamlike  Function,  \ , , Contours. 
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Streamlike  function,  x^>  Contours 
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Fig.  3.  Secondary  Flow  Velocity  vectors 


Fin.  7.  Streamlike  Function,  x 
y on  the  Upper  Side. 
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Fig.  8.  Streamlike  Function,  x 
v on  the  Lower  Side. 
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Streamlike  Function,  xi'  Contours 
x on  the  Right  Hand  Side. 
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